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The tunnel magnetoresistance (TMR) of F/O/F magnetic junctions, (F's are 
ferromagnetic layers and O is an oxide spacer) in the presence of magnetic impu- 
rities within the barrier, is investigated. We assume that magnetic couplings exist 
both between the spin of impurity and the bulk magnetization of the neighboring 
magnetic electrode, and between the spin of impurity and the spin of tunneling 
electron. Consequently, the resonance levels of the system formed by a tunneling 
electron and a paramagnetic impurity with spin S = 1, are a sextet. As a result the 
resonant tunneling depends on the direction of the tunneling electron spin. At low 
temperatures and zero bias voltage the TMR of the considered system may be larger 
than TMR of the same structure without paramagnetic impurities. It is calculated 
that an increase in temperature leads to a decrease in the TMR amplitude due to 
excitation of spin- flip processes resulting in mixing of spin up and down channels. 
It is also shown that asymmetry in the location of the impurities within the barrier 
can lead to asymmetry in I(V) characteristics of impurity assisted current and two 
mechanisms responsible for the origin of this effect are established. The first one 
is due to the excitation of spin-flip processes at low voltages and the second one 
arises from the shift of resonant levels inside the insulator layer under high applied 
voltages. 

1. Introduction 

The observation of the large tunneling magnetoresistance effect at room temperature in 
tunnel junctions of the form M/O/M' (where M and M' are magnetic metals and O is 
an oxide tunnel barrier) has stimulated a renewed interest for these systems 0, [§. 
Besides the fundamental interest for spin-polarized transport, these structures are also 
foreseen as potential candidates for sensitive magnetic sensors and memory cells in random 
access memory devices. The first model of spin-dependent tunneling in the framework of 
classical quantum mechanics was proposed by Slonczewski |4j] . However, in this approach 
no scattering of electrons in the magnetic metallic electrodes was taken into account. This 
model has been subsequently developed in Refs. f|, H by using the Kubo formalism of 
linear response. The effects of elastic impurity scattering inside the metallic layers and 
at interfaces between the dielectric and conductive layers could then be incorporated in 
the model. On the other hand, it is well known [0] that the presence of impurities inside 
the potential barrier can lead to the mechanism of resonant tunneling when the localized 
electronic states within the gap of the insulator formed by embedded atoms lie close to the 
chemical potential of the system. This situation was qualitatively studied in mesoscopic 
semiconductor system || in case of one- and two-impurity resonant channels by means of 
classical quantum mechanical treatment. The same approach afterwards has been used 
in [[J with the application to the impurity-assisted tunneling magnetoresistance (TMR). 
The numerical analysis of this problem which was carried out in Refs. [pi], |TTJ should 
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also be mentioned. In paper |9j only the case of spinless impurities was considered, 
and the author came to the conclusion that the TMR amplitude decreased due to the 
impurity assisted tunneling. The problem of the paramagnetic impurity assisted tunneling 



in tunnel magnetic junctions was investigated recently in Ref. ||12|| . The author claim, 
that he investigated resonant tunneling through impurity resonance level, however he did 
not introduce the line-width of this resonance, which, as it will be shown below, does 
depend on the position of the impurity atom inside the barrier as well on the magnetic 
configuration of the magnetic layers. As it will be shown below, namely these line-widths 
define the value of the tunneling conductance and the amplitude of the TMR for spin- 
conserving and spin-flip resonant tunneling. An attempt of the analysis the same problem 
has also been undertaken in Ref. [y|, but nevertheless the microscopic mechanism of 
electron scattering on the paramagnetic impurity has not been taken into account. 

In this paper, we propose a renewed study of the problem of impurity-assisted tunnel- 
ing in spin- valve junctions of the form F1/O/F3, where F's are ferromagnetic electrodes 
and O is an insulating barrier with embedded paramagnetic impurities, that incorporates 
the effect of both elastic and non-elastic spin-flip scattering due to the exchange interac- 
tion between the itinerant electrons forming the tunneling current and the localized spins 
of impurities. It will be shown that non-elastic scattering has an essential impact not 
only on the temperature variation of the TMR (which is a well established result |1|) but 
also on the I-V characteristics of considered structures. The latter effect was predicted in 



Ref. [IS], where the TMR dependence on the electron scattering on interfacial magnons 
was investigated. 



2. Model 



2.1. Kubo formula and general expression for the conductivity of the system 

The following simplified model is adopted throughout the paper. First of all, the thickness 
of an oxide layer is supposed to be much smaller than its in-plane dimension, so that the 
system may be considered as homogeneous in the xy-plane (parallel to the interfaces) 
and inhomogeneous only in the ^-direction (growth direction). Within each layer, the 
electrons are described as a free-electron gas and they undergo scattering on the 3-D 5- 
function impurity potential within the insulating barrier. Within these approximations, 
the Hamiltonian of the system has the form 

H = H + H int , 

where 

Ho = -^y A + U ( z ) ~ ^ B Hf(z)(s z + S z ) (1) 

flint = a 6 ( T ~ C i £ - J ( SS )I • 

i 

Here the summation is performed over the location of impurities Cj inside the barrier, 
ao is the lattice constant, e denotes the scattering potential amplitude on the impurity, 
J is responsible for the s-d type exchange interaction between a conduction electron spin 
s and the impurity spin S, U(z) is a model step-like potential seen by the conduction 
electron as it is represented in Fig. 1. We take into account the exchange splitting of 
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Fig. 1: The potential profile seen by electron propagating through the F/O/F junction com- 
prising impurity defect inside the oxide spacer, ki, k%, q2 are the momenta inside the magnetic 
layers and oxide barrier, respectively. V^j^ denotes the spin-dependent conduction band bot- 
tom, U is a level of the barrier and ef is Fermi energy. The paramagnetic impurity is located at 
point c. The variation of the potential profile under high bias voltage is indicated by the dashed 
line. 

the (i-band by introducing different values V^ 1,3 for the position of the bottom of the 
conduction band in Fi and F3, depending on the mutual orientation of magnetization 
in the layers and the spin \i =|, j of the conduction electron. Hl s (z) represents the 
effective field acting on impurity and electron spins inside the barrier. The origin of this 
field is the super-exchange between the spins in the bulk of ferromagnetic layer and in 
insulating layer. We suppose that H\ (z) decreases exponentially with the distance from 
the interface in the depth of the oxide layer. m(z) corresponds to the effective electron 
mass that we suppose is equal to m in the ferromagnetic layers and to mo in the insulator. 
Throughout the paper, it is expressed in units of bare electron mass m e . We also assume 
that the mass of free-like electrons in ferromagnet is slightly differs from m e , i.e. m « 1 
and we will eliminate it from all subsequent expressions. 

We start from the Keldysh technique for Green functions together with Kubo exact 
formula of linear response theory for the static conductivity which relates its real part 



with the current-current correlation function and may be written in the form |Lc 

r+00 



2k B T 



j p {v\t')j^v,t))d{t-t') 



1 f ^L) l^( G i P (^y^)^rVr'G>^r\t',v,t))d(t-t'), (2) 



2k B T \2m) 

<— > — > < — 

where /1, p denote the projections of the spin of the electrons, V r = ^(^r — ^ s ^ ne 
asymmetric gradient operator and G< p and G> are corresponding Green functions in 



Keldysh formalism [17j]- (• • •) represents the quantum statistical averaging over the dis- 
tribution of impurities and degrees of freedom of impurity spin. This expression is most 
general and holds both for the elastic impurity and defect scattering or inelastic, includ- 
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ing magnon and phonon, scattering. To evaluate expression (|2|) one needs to introduce 
the retarded Green function G„„( ) that in our case is defined by the following 

differential equation: 

f + £1? - ^ - u(z) ] G >' z '' K - £) = s " S{z ~ z ' ] 

in the mixed real-space momentum representation |5|, ||, where k = (k x , k v ) is the com- 
ponent of the electron momentum in xy-p\a.ne of the layers and z is the coordinate per- 
pendicular to xy-plane. We should note, that by definition the conductivity (|2|) is defined 
as a linear response on the externally applied electric field and does not depend on z and 

Oj(z^) 

z' because of the obvious condition — — = 0. 

oz 



Let now denote k% — J 2(e — Vi), k% = J 2(e — V3 1 ) momenta of electrons with energy 

e and spin /x in the ferromagnetic layers and g = ^2m (U — e) is an imaginary momentum 
inside the barrier. By introducing the following functions 0111 = ft/go'- 



d£(x) = \jk^ 2 - g$x 2 , c$(x) = \fk% 2 - g 2 x 2 , g 2 (x) = g Vl + x 2 , 

the final expression for the conductance of the system, comprising only one impurity, 
located at point c, at given temperature T is written as follows: 

a(T,c) = a (T) + a imp (T,c). 

The first term is given by 

g 2 e 2 [+°° f df(e)\ f< xdx W^mtgle- 2 ^ 
ao[1) toh^J-co [ de J Jo 2n (m 2 cf + g 2 )(m 2 cf + g 2 Y [ ) 

where Xq = min |^-, ^ j, f(e) = 1 + e^ £ ~ £F ^ is Fermi function, and w = b — a is the 
width of the insulating spacer. It represents the well known result for the pure tunneling 
conductance ||. The second term cx imp (T, c) is directly related to the impurity assisted 
tunneling. It is convenient to write it down as a sum of two contributions: 

where the first term corresponds to the conductivity due to elastic spin-conserving pro- 
cesses of scattering electron on the impurity and the second one summarizes all other 
events when the conduction electron changes its spin after tunneling through the barrier. 
We have derived the analytical expressions for these two terms, which are valid under two 
assumptions. Namely, under domination of single electron scattering on impurities over 
multiple scattering of two and more electrons on the same center and under absence of 
polarization of impurity spin induced by the ejection of spin-polarized electrons. Then the 
final result for these terms is written as (the details of its derivation are outlined further) 

ff?P(r ' c) = S (If) C & ~dT" f) ((g( £ » tf »K(c)*f (e)- 

dfi(e + fi B Hf) ^mf^KLf^R, 
de 



(ti(e)yti(e))^(c)^(c) , (4) 
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<T P (T, c) = - ^LJ — J_2 de{Me - fi B Hf) [l - h{e + p B Hf)] (L(e)t + (e)) + 
fl(e + nsHf) [l - Me - VsHf)} (f + (e)L(e))} x i {$f (c)$f (c) + $f (c)$f (c)} . 

Here S is the junction area, $^^(c) and $^(c) are the probabilities of tunneling of the 
electron from the left or from the right electrode to impurity, located at point c. Omitting 
the exponentially small terms, the expression for these probabilities can be written as 



L xdx 2djmlql 

" ' 2vr Kef + ql) e 



-2q 2 {c-a) 



*>R M = f <BX xdx 2c WA -2g 2 (6-c) 

Pl J " A) 27T Kef + g2) 

The quantities ((tl^)WJ-^)(e) and (£_£+) (e), (£+£_) (e) in (£|) represent the scattering 
amplitudes of electron on the impurity center for the case of spin conserving (\in, |) — > 
\out, t)) or |m, |) — > | out, j)) and spin-flip (|m, j) — > \out, j) or |m, j) — >■ |owt, |)) transi- 
tions averaged over the distribution of paramagnetic impurity spin. Here \in) and \out) 
denote the initial and final states of impurity. Operators and t± form a one-center ma- 
il 

h a 

spins and are expressed as 



trix t = ( j z *h ) in the direct product of the linear subspaces of electron's and impurity's 



tl (i) (e) = ™i VP(s), 

£ / \ \ g a\J jl . . 

± \-V^\e)G m {e) ± l-a 3 (e + lJS z )G n) (ey 

where effective potentials Vj^ are given by 

VP(e) = a 3 L T -JS Z + -S T "MM^ § d 

Here Gm)(e) is the electron Green function at point c: 

/" Kmax ndn 
G M (e,c) = / G^ K (e,c) ——, 

JO Z7T 

where K max = is a cut-off of in plane momentum that stems from the finite size of 
Brillouin zone (we substitute the Brillouin zone's projection onto (k x , k v ) plane by the 
circle of radius K max of the same square in Ky-plane). The real and imaginary part of 
G%(e, c) (/i is the spin index) in the leading order of magnitude are given by 

ReG^e,c) = -^, (6) 
1i 



Im^( £ ,c) = -^(c)+$«(c 

Let us now explain the derivation of expression (H) and clarify the two assumptions 
under which this formula is valid. To derive (f|) from the starting point (@) one can first 
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Fig. 2: Two diagrams that make contribution to the conductivity in the second order of per- 
turbation theory in case of a) elastic spin- conserving scattering and b) non-elastic spin-flip 
scattering. 

of all examine two diagrams (a) and (b) (see Fig. 2) that contribute to spin-conserving 
and spin-flip part of <r imp (T, c) at second order of J, respectively. One may easy verify 
that a general structure of these diagrams is just the same as the final result in form (||) 
with the mere difference that the one-center t-matrix is reduced at first order of J to the 
initial potential Eq — y f* ^ • 

Moreover, the diagram (b) contains both direct and indirect processes in equal pro- 
portion with common factor 1/2 for any of the possible channels |in, f) — > \out, I) or 
| in, |) — > \out, |). The thermodynamic averaging (. . .) in the second-order expansion is 
simply reduced to the averaging over the Boltzman distribution of the impurity spin in 
the "external" effective magnetic field H° s which was introduced in ([!]), i.e. with the 

density matrix p = Z^ 1 exp j 2ilB ^ T ^ z j, where Z = 2 cosh ( 2/ ^y z ) + 1. After that, 

it is easy to check that the total probabilities (with account of Fermi factors of electron 
states) of direct and inverse processes are equal which means that the principle of detailed 
equilibrium holds. In particular it leads to the vanishing of spin-flip processes in a system 
at zero temperature and vanishing voltage bias. 

After that preliminary discussion two assumptions should be made to justify the re- 
sult (§: 

i) We assume that the occupation of the given impurity center simultaneously by 
two electrons with different spin (due to Pauli principle) is a rather rare event or, in 
other words, we do not take into account many-electrons effects. It may be justified: a) 
by Coulomb interaction between electrons that make unprofitable their arrangement at 
the same site of the lattice; b) by the large number of impurity centers that provides 
a sufficiently large number of one-step channels so that electrons may be considered as 
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independent. 

ii) We also neglect the influence of electron current on the statistical distribution 
of paramagnetic spins inside the oxide barrier. This assumption is valid for practical 
intensity of tunneling current which is low enough not to produce a spin polarization of 
impurities by injection of spin-polarized charge carriers. 

Under these assumptions the expression can be obtained by simple substitution of 
scattering potential H mt at site c, on Fig. 2 by the corresponding one-center t-matrix in 
accordance with (|5|) and assuming that the averaging over the degrees of freedoms of the 

impurity is carried out by means of unperturbed density matrix p = Z" 1 exp j ^^Jr ^ 

In this form the structure of the result (^) is similar to the one obtained in Ref . [pl| , where 
the spin-flip scattering of electrons at interfaces of tunnel junctions was investigated in 
the framework of tunneling Hamiltonian and the second order perturbation theory. 

2.2. Resonant tunneling in the case of nonmagnetic impurities 

To extract the physical nature of resonant tunneling through the impurity states contained 
in expression (^) we proceed as follows. For the sake of clarity and simplicity we consider 
first the case of zero-spin impurity. Then only one element of t-matrix at cite c survives 

a3 ° £ ° 



1 - age G T (;)(e)' 

It defines the position of a resonant level inside the gap of the dielectric band structure 
by finding the root of the equation ale ReGm)(ei) = 1. From expression @, it follows 
that the real part of the Green function Re G^e, c) is independent on c and spin /z up to 
exponentially small terms. Therefore, the position of level Si is weakly depend both on 
the position of impurity inside the barrier and on the direction of the spin of tunneling 
electron. Evidently, only those impurities for which Si is close to the chemical potential 
Ep contribute to a significant extent of the total current at low bias voltage. Therefore, it 
is possible to expand the denominator in t (e) in powers of (e^ — e) . If we now introduce 
the position dependent line-widths 

r£(c) = $;(c)/ReG'M, rJ(c) = $J(c)/ReG'M, (7) 



where ReG'(e F ) = £-ReG{e) is the energy derivative of the electron Green function 



at Fermi level, then we obtain the general formula for the resonant case of impurity 
assisted tunneling 

where r^(c) = T^(c) + r^(c) and the summation by i is performed over all resonant 
levels. For the qualitative analysis, one may evaluate the expressions (0) for T^ L )(c) 
approximately by considering the case k = which is valid if e~ 2qoW C 1. In this 
approximation 
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and expression (||) reproduces the result of Ref. ||. To proceed further, we discuss some 
assumptions concerning the parameters of the model. We consider the case of Co electrode 
and AI2O3 as an oxide layer and take typical values of k^ = 1.09 A -1 , kf = 0.42 A^ 1 , 
m ~ 1 for itinerant electrons in Co and a typical barrier height for AI2O3 (measured from 
the Fermi level Bp) U — Sp = 3eV with an effective mass m = 0.4 (Ref. ||), that gives 
q ~ 0.56 A -1 . Assuming the thickness of the barrier w ~ 20 A -1 , one may estimate the 
conductance o"o of the system without impurity from (|3]) by means of approximate formula 

°-**h\w)% (m*k% + qt){mlk™ + ql) ' 1 UJ 

that leads to GMR ~ 16 %. To estimate the value of the line-width (^) we consider 
impurities located close to left interface at a distance, say, of two atomic layers which 
corresponds to (c — a) ~ 4 A. For spin up electrons it gives r-j-(c) ~ 0.02 eV. Further 
in this paper we restrict ourselves to the case of temperature interval from 4.2 — 300 K 
(0.025 eV). We assume that the impurity levels Si in the band gap of the insulator form 
a narrow impurity band of width Ae which spreads symmetrically with respect to Fermi 
level Ef and, following Ref. ||, we introduce its density of states v(e) per unit volume 
and unit energy interval. We assume that Ae is of the order 0.1 to 0.2 eV, i.e. an order 
of magnitude greater that the above estimated line- width. In this context, with a good 
accuracy, the impurity conductance (H) rewrites as follows: 



a imp < 



(c) - T u(e P ) £ \- Ts j p(e, c) de, (11) 

where factor 



, s 2 f Ae 

p[e, c) = — arctan 



V 2r V( C )/ 

arises from the integration of exp. (]8|) over impurity levels in the range of impurity 
band. Due to above mentioned estimations equality p(e, c) ~ 1 holds with a good degree 
of accuracy and in this case exp. ( |TTD becomes in agreement with Refs. [§. 0. 

2.3. Resonant tunneling in the case of paramagnetic impurities 

To investigate the general case of paramagnetic impurity we follow the same procedure 
as in the previous section. Let J = s + S be the total magnetic moment of the system. 
We may state that [H, J z ] = and, therefore, J z is a good quantum number. We regard 
the t-matrix (|5|) as an operator acting on the spinor subspace \o~, m), where a = ±~ and 
m = ±1,0 corresponds to the projection of the ^-component of electron and impurity 
spin, respectively (we consider the case S = 1). As far as its total magnetic moment 
along the z-axis J z = s z + S z is conserved, the matrix elements (o~imi\t\o~2Tri2) are non- 
zero only if mi + o~± — ni2 + 02 and, therefore, it is convenient to introduce the notation 
tm° 2 = ( cr i m iN cr 2^2), where rrij = mi + o\ = m 2 + 02- These elements are simply 
calculated from (|5|). The non-zero ones are written as follows 

a 3 (e - J/2) 



t 3/2 



l-o§(e - J/2)G t (e) 

al(e - J/2) 
l-al(e Q -J/2)G i (ey 



AX _ "-pycu - 01*) / 10 x 
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and £±i/2 = ( fj 1/a 4 f T 1/2 J corresponding to the subspace rrij = ±| with 



"±1/2 "±1/2 



*i) 2 (*-i/ 2 ) = -^^{eo + J/2-a 3 G Hl) (e)(e -J/2)(e + J)} 

^ 2 (t U i/ 2 ) = 7r ^{eo-agGi(T)(e)(eo-J/2)(eo + ^)}, 
^±i/ 2 (£j k J 



(13) 



t 



J, 



^ ±1/2 " ±1/2 " v/2A ±1/2 
where denominators are 



A 



±1/2 



£ 



(l - agG i (e)(e - J/2)) (l - aj>Gr t (e)(e + J)) ± J(G T (e) - G^e)). 



As can be seen, two poles of the t-matrix defined from equations a[]Re £7(63/2) (£0 — J/2) = 
1 and aQReG(ei/2)(£o + J) = 1 correspond to two multiplets £3/2 and £1/2 with a total 
angular momentum j = 3/2 and j = 1/2, respectively. If J > 0, then e 3 / 2 < £1/2, i-e. the 
multiplet with j = 3/2 has a lower energy than one with j = 1/2. As for non-magnetic 
impurity, we restrict ourselves by considering the regime of only one-channel resonant 
tunneling. We assume that J > and the lowest impurity levels £, = £ 3 / 2 corresponding 
to the multiplet with j = 3/2 lie close to Ep . We note that the typical value of exchange 
coupling J is of order 1 eV and due to this fact we may eliminate the resonant level 
£1/2 from further consideration. Then, as in the previous analysis for a non- magnetic 
impurity, only the resonant part of the t- matrix (|i~2"l , |T3"l) at energies close to chosen £3/2 
is essential for the subsequent calculations. Expressions (|12]JT3|) can be easily written as 
follows: 



C ±3/ 2 \ E ) 
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(14) 



where 7j(c) = §r-f(c) + |Tj(c), 7;(c) = 3^(0) + |Tj(c) are the inverse lifetimes of the 
resonant states with rrij = ±1/2. This result allows simple qualitative interpretation. Let 
us look, for a example, at quantum states with rrij = 1/2. From elementary quantum 
mechanical theory one may conclude that 



l) l/2 z 
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: 2' m i 
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(15) 



As we have assumed, only 
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2> m j 



gets into resonance and, therefore, e.g. 



t 



n 
1/2 



'1/2 



^i/2 



y/2/3 in agreement with (0). On the other hand, 
T,m s = oWa/I l,m s = 



2 ,rrij — 2 
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and, hence, its inverse life time is given by t x i 2 = 1 + |tj . 

We substitute all t-matrix elements in (|) by its resonance expansion (|r3D. To proceed 
further, one has to perform in (J|) the configuration averaging over all impurity centers 
and thermodynamic one over all possible channels. Suppose, that the impurities are 
distributed uniformly in the space in the interval [zq — A, zq + A] along z-direction with 
the width of 2A and center zq which we have chosen to be close to left (L) ferromagnetic 
contact. After that, it is possible, first of all, to average the Lorentzian peaks over the 
distribution of impurity levels by averaging them over £ 4 in the range of impurity band 
with factor v{£f) and to perform the thermodynamic averaging by integrating over e 
and neglecting the dependencies of T L ^(c) on energy. On the third step, the averaging 
over the space distribution of impurities along z-direction should be made. Following 
the outlined procedure, the total conductance @ as a function of temperature for the 
parallel and antiparallel alignment of magnetization of the ferromagnetic layers is written 
as a sum of factorized terms over all possible scattering channels: 

fipmj 

a A P{T) = 2 -fY.P% (^F) AP < P ^o, A)»(e F ) + of, (16) 

fiprrij 

where 

1 rz +A 

<J (z , A) = — y^_ A o% (c)p(e F , c) dc. 

The origine of p(ep,c) is the same as in exp. ( pTTj) and explicit form of functions P^f.(h) 
and er^.(c) is given in Appendix A. We have also used the same notation of matrix 
indexes as it was previously introduced for t-matrix elements. o~ P and (Tq P are the tunnel 
conductances of the pure system in accordance with (Q). Factors Pfif. and a^. represent the 
thermodynamic and quantum mechanical probabilities of the given process, respectively. 
Expressions (|IE|) are the final results of this section and their analysis is presented below 
(see Sec. III). 



2-4- Dependence of conductivity on bias-voltage 

We are also interested in I(V) characteristic of the considered system. To derive the 
general formula for the current, one may simply extend the expressions (f|) to the case of 
finite applied bias voltage. Consider, for example, the contribution to the total current 
J, coming from all possible channels of the form \in, f ) — > \out, j) for tunnel electrons 
moving from the left electrode to the right electrode and of the form \in, j) — - > \out, t) 
for electrons moving from the right to the left, respectively, i.e. in both cases an electron 
has an "up" projection of spin in the left contact and "down" projection of spin in the 
right one after or before scattering. From the general concept, one may conclude that 
this contribution to the current can be written as 



I n (V) 



1 2e 




S \ nh 



de\fi(e-fi B Hf-eV) l-f^e + ^H] 



TCS\ 



t(e)U(e))- 



f^e + psHf) 1-Me-» B H? - eV) (i + (e)L(e)) W (c)$f (c),(17) 
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where it is assumed that the voltage bias is applied from the left to the right direction. 
It is important to notice that inelastic spin-flip processes of the electron scattering on the 
impurity were taken into account in derivation of the exp. (O) but they were omitted 



in Ref. |12| . Analogous expressions can be written for all other channels. In the case 
under consideration, expression (|TT| ) contains two regimes of non-linear behavior of I{V) 
characteristic. The first one reproduces a zero bias anomaly due to excitation of spin- 
flip processes at low bias voltages of order of magnitude /ieif| ff (we believe that it is of 
order 5 mV). In this range, as before, one may assume that the resonance amplitudes 
{(tl^Y (tl^)) and (t-t + ) are nearly independent of the energy after averaging over all 
possible configurations of impurities. As a result, the voltage dependence of total currents 



for parallel and antiparallel configurations are given by formulae similar to ([Tq) : 



I P (V,T) = | £ I%(V,Hf) p a%(z ,A)v(e F )+ a p V, 

lipm-j 

I AP (V,T) = ^J2l%(V,Hf) Ap a%(z ,A)v(e F )+a£ p V. (18) 

The expressions for I^XYi Hf 1 ) are given in Appendix B. The voltage dependent conduc- 
tances (J P (V, T) and a AP (V, T) can be obtained from ([17D by derivation with respect to 
V. The detailed analysis of this physical situation is presented in the next section. 

The second source of possible non-linear character of I(V) dependence is the variation 
of potential profile U(z) (see Fig. 1) under applied bias voltage. It follows, then, that 
the latter introduces the correction to exp. (|9]) and ([H]) and they can be calculated 
with the use of Wentzel-Kramers-Brillouin (WKB) approximation jTj| assuming that the 
applied voltage produces the uniform electrical field inside the insulating layer. In the 



case of pure tunnel conductance it is known[19] that both conductances for the parallel 



and antiparallel configurations increase with the increase of applied voltage so that the 

IP(V) - I AP (V) 

TMR as a function of V, defined as ttftttn > drops significantly at the voltages 

I Ay (y) 

of order 1 eV. The contribution of impurity assisted tunneling may change considerably 
this situation in the case of non-uniform spatial distribution of impurities, e.g. when they 
are distributed in the vicinity of only one electrode. In this particular situation, as we 
will show, the essential variation of TMR amplitude in the case of magnetic impurities 
(in contrast to non-magnetic ones) will take place at the bias voltages compared with 
impurity band width Ae. 

For the sake of simplicity we consider, first, the case of non-magnetic impurities. In 
the WKB approximation the contribution from all impurities, located at given point c, 
to the total current I(V) will have the form similar to (§) and fllTl) : 

3 

where 



imp (c) = T v(e F ) £ /_ m {f(s - eV) - /( £ )} -g^^ p ( e , V) de, (19) 



r L( r \ = k l^a m T a 1 -Sg/h. 
MU (9a) 2 +Af>f ' 

r?(c) = ( k kT m °S 2 ^ /n - (20) 
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Here ql = q$ = 2mo(U — e), q 2 a = q 2 . + 2moeV are imaginary momenta of electron with 
the energy e in the vicinity of the right and the left electrode, q^ = qo ± - — ^ — ~j E is 



the electric field in the device. We also introduce q c = q + 2m eV(b — c)/w imaginary 

3 3 3 3 

momenta of electron on the impurity center. Then S a = — — , S& = — — represent 

3moeE 3moeE 

the classical actions along the path from the left contact to the point c in the barrier 

eE 



and, afterwards, from this point to the right contact, respectively; r a = ^ a and 



n = ^ c ^ b denote the passage times associated with these paths. Factor 
eE 



p(e, V) — — { arctan 

7T 



As 
2 



r» + r*(c) 



arctan 



e F 



2 



r^(c) + r*(c) 



as before, arises from the summation over all impurity levels £j and gives the relative 
weight of all resonant channels with energy e. To clarify the situation, it is sufficiently 
to consider the most resonant channel with energy e T = Bp + eV(b — c)/w at which 
p(e, V) reaches its maximum. One may note that e r corresponds to the resonant impurity 
level that exactly coincides with Fermi energy at vanishing voltage and it shifts linearly 
with the increase of applied bias depending on the position c of impurity inside the 
barrier. As it was stated earlier, the most interesting case takes place when the point c 
is situated close to the left contact. Then one can see that T L (c) > T R (c) and, thus, 
j imp (c) ~ T R (c)p(e r , V). At bias voltages much more lower than the height of the barrier 
ip = (U — £p), Sb can be expanded in powers of V: 

which shows that r iJ (c) ~ exp(— Sb/fr) is an increasing function of V in the vicinity of 
V = 0. Hence, it leads to increase of differential conductivity cr(V) = dl/dV under direct 
bias voltage, applied to the barrier from the left to the right direction, and to decrease 
of cr(V) under inverse bias voltage. The physical meaning of such behavior is rather 
obvious. From expression for S& it follows that electrons tunneling under forward bias 
due to resonant levels lying close to e r will propagate through the potential barrier with 
the height which is effectively less compared with that in case of inverse bias. 

The expression for the paramagnetic impurity assisted current at finite voltages has 
the structure similar to exp. with Fermi distribution factors written in accordance 
with the general formula ( |T7D and integrand expression has the form given in Appendix A, 
where line- widths I^q have to be substituted by WKB approximation (p0|). In the case 
of magnetic impurities the above outlined mechanism of asymmetry in I(V) characteristic 
due to the resonance levels will essentially contribute to the voltage dependence of TMR 
in question, and it will be discussed in the next section. 

3. Results and discussion 

In this section we consider the temperature and bias voltage dependencies of the con- 
ductances and TMR effect of the considered structures. We investigate the case of 
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Fig. 3: Tunnel magnetoresistance at T = as a function of polarization P = — p p- under 

fixed A;^ = 1.09 A^ 1 . Other parameters are: go = 0.56 A" 1 , mo = 0.4. Solid line corresponds 
to the case of impurity concentration x = 8 x 10~ 5 , dashed line represents the case of absence 
of impurities. 

Co/Al 2 3 /Co junction with the typical parameters that have been already mentioned 
in section 2: k^ = 1.09 A -1 , k^ = 0.42 A^ 1 are the Fermi momenta of itinerant electrons 
in Co, go — 0.56 A -1 is the imaginary momentum in the barrier, mo = 0.4 is the effective 
mass in the insulator and w = 20 A is its thickness. We focus on the most interesting 
situation when impurities are introduced at the vicinity of the left electrode at a depth 
w\ inside the insulator layer. We chose width w\ = 4.06 A that corresponds to two atomic 
monolayers. The essential parameter of the model that has to be defined is the effective 
molecular field if| ff acting on impurity spins. One may suppose that it should exponen- 
tially decay in the depth of the barrier. We have, therefore, set it to HbH% = 5 meV 
(58 K) that is of two order less than the critical temperature in the bulk Co. 

Consider, first, the case T = 0. It is possible to estimate the concentration of impu- 
rity atoms so that its contribution to the resonance conductivity is comparable with the 
ordinary tunnel conductance of the system. One can write that the impurity density of 
states (see exp.(|Pl|)) v^Ef) = g^/^g) ' wnere ^ * s a total number of impurities and Ae 

is the width of its energy distribution. On the other hand N{ = xN and iV = Swi/a^ 
where N is the total number of atoms in the layer which contains the impurities, x is 
the local concentration of impurities in this volume and ao is the lattice constant. This 
yields v(ef) = x ^ ■ We introduce the characteristic concentration xq defined so 

that in the case of nonmagnetic impurity, the impurity conductance (0) is equal to the 
tunnel conductance ( P"U|) of spin f channel in the parallel magnetic configuration of the 
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ferromagnetic layers. Such definition leads to 

_Ae/ 2m ag \/ w \ / a \ kf q m (q oWl ) exp (-2g Wi) 
^"vr 2 ^ ft 2 J U - w J \w J kf 2 ml + ql 1 - exp (-2 9o Wi) 1 J 

If we choose Ae = 0.2 eV, then xq = 6.5 x 1CT 5 . The conductance of the system at T = 

can be extract from the general expression (14). We suppose that for both parallel and 

antiparallel configurations the left electrode has "up" magnetization and, hence, H%r is 

positive in both cases. At zero temperature, all spin-flip processes are frozen and due to 

the above assumption, only the configuration of impurity spin with m s = 1 is possible. As 

a result, only two resonance channels from many possible ones have nonzero contribution 

to conductivity, namely | ] ) m s = 1) «-> | ] 1 m s = 1) and | j, m a = 1) «-> | I, m s — 

1) with rrij = 3/2 and 1/2 respectively. From exp. (Al) (see Appendix A) it follows 

that the channel with rrij = 3/2 gives the main contribution into the conductivity at 

low temperatures and P (jlj 2 ~ ^ an d AP<7 3/2 ~ ^l- So these contributions depend on 

the mutual orientation of the magnetization of the ferromagnetic layers, therefore they 

p _ AP 

increase the total amplitude of the TMR. The total expression for TMR = — — -^p — 
including all possible channels may be written as 

R = r r r J (2 + i(i + ir,/7,)) (22) 

where x is concentration, xq is defined by (pl|), 7-f = |r| + |r^ and = k^q^ / (k^ml+ 
ql) are the tunneling density of states for j (J.) spin electrons. The dependence of TMR 
effect versus the polarization P = (k? — kf) / '(kf + kf) is shown in Fig. 3 in comparison 
with the non- resonant tunnel conductance at x = 0. The total TMR amplitude is larger 
than the TMR due to direct tunneling, in accordance with considerations written above. 
For a given polarization P = 0.44 in case of chosen parameters, the contribution of the 
impurity assisted tunneling leads to strong enhancement of TMR amplitude (typically by 
a factor 2, see Fig. 3). 

We note, that in the case of nonmagnetic impurities, distributed in the vicinity of only 
one contact, the resonant impurity conductance <r imp ~ + is equal for both parallel 
and antiparallel configurations and, therefore, in this case the mechanism of impurity 
assisted tunneling is not able to enhance the TMR effect. The enhancement of TMR 
amplitude in the case of paramagnetic impurities is essentially due to the presence of fer- 
romagnetic exchange coupling between the magnetization in ferromagnetic electrode and 
the impurity spins which tends to induce a ferromagnetic order in the plane of impurities 
and, as the result, leads to the preference of impurity spin to be found in the quantum 
state with m s — 1. 

The temperature dependences of resonant conductances for parallel and antiparallel 
configurations in the interval from 4.2 K - 300 K are presented in Fig. 4. In the case of 
parallel alignment, a^iT) is nearly independent on the temperature, but in antiparallel 
situation there is a 50% increase of impurity conductance a^f (T). This originates from 
the thermal excitation of both spin-flip and spin-conserving processes which are frozen 
at zero temperature. For AP configuration the process | ].,m s = 0) —>■ \ i,m s = 1) was 
forbidden at 0°K but now it is allowed and gives a large contribution into the current as 
it is proportional to the product of the largest density of states IYTf. As a consequence, 
the TMR effect decreases with the increase of the temperature. 



14 




Fig. 4: Tunnel magnetoresistance as a function of the temperature. (On the insertion: con- 
ductance a for the parallel and antiparallel configurations.) Parameters are: k^ = 1.09 A -1 , 
kf = 0.42 A -1 , qo = 0.56 A" 1 , tuq = 0.56. Concentration of impurities x = 8 x 10~ 5 . 



We have also calculated the dependence of the differential conductances on bias voltage 



according to (|18|) and (|19l) . These dependences at T = 4.2 K and T = 77 K are presented 
in Fig. 5. A new effect is predicted: the voltage dependence of the conductance in the 
antiparallel alignment of magnetization in ferromagnetic layers is asymmetric under for- 
ward and inverse bias voltage when the paramagnetic impurities inside the insulator layer 
are distributed close to only one of the interfaces and are bound by exchange interaction 
with magnetization of the nearest ferromagnetic layer. 

One can distinguish two different mechanisms that give rise to the presented asym- 
metrical behavior with respect to inversion of bias voltage. The first one, which we refer 
as zero bias anomaly, manifests itself at low voltages of the order of 10 mV (for particular 
chosen parameters in our model) and is strongly pronounced only at low temperatures 
(see Fig. 5, the case of T = 4.2 K). It originates from the excitations of spin-flip processes 
on the impurity centers. One may look at the general expression and consider the 
case of low temperature. An electron undergoing spin-flip scattering, may transfer an 
amount of energy uj = l^i^H^ to the impurity spin thus exiting it at an higher energy 
level or on the contrary may acquire this quantum of energy from it. The latter process is 
impossible at low temperature. The former one is possible only if an electron moving, say, 
from the left contact possesses an excess energy of at least ioq with respect to Fermi level 
in the right contact. The only one process that contributes to this anomaly at low tem- 
perature is 4>\i 2 — > (see exp. [TJ). For antiparallel alignment of the magnetization its 

quantum mechanical probability is proportional to AP 0"jy 2 ~ frJ/TT for electrons moving 

from the left ferromagnetic layer into the right one and is proportional to AP <7y 2 ~ 1^/7? 
in the case of electrons moving from the right to the left. For parallel configuration of 
magnetizations these probabilities are equal in both directions and are proportional to 
P<T i/2 = P<T I/2 ~ §^1-^1/71- A s a resUi t, zero bias anomaly at T = 4.2 K looks asymmet- 
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Fig. 5: Differential conductance as a function of the bias voltage at T = 4.2 K and 77 K. 
Parameters are: kf = 1.09 A -1 , kf = 0.42 A -1 , q = 0.56 A -1 , m = 0.4. Thick dashed 
and solid lines correspond to the conductance in the parallel and antiparallel configurations 
respectively, calculated in the WKB approximation. For the comparison, the same dependences 
at T = 77 K, calculated by approximate formulae, are indicated by thin lines. Concentration of 
impurities x = 8 x 10~ 5 . 
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Fig. 6: Tunnel magnetoresistance as a function of the bias voltage: solid line - T = 4.2 K, 
dashed line - T = 77 K. Parameters are: kf = 1.09 A -1 , kf = 0.42 A -1 , q = 0.56 A -1 , 
mo = 0.4. Concentration of impurities x = 8 x 10~ 5 . 

rically in the case of antiparralel configuration and is symmetrical in the case of parallel 
alignment of magnetizations. 

The differential conductances as a functions of the bias voltage at T = 77 K have been 
calculated using two different approximations. Thick dashed and solid lines correspond 
to the conductances in the parallel and antiparallel configurations, respectively, that have 
been calculated by means of WKB approximation in accordance with expressions ( |l7|) 
and fll9p . For the sake of comparison, the same dependences, indicated by thin lines, 
have been calculated with the use approximate formulae fll8D , where the dependence of 
t-matrix elements on the applied voltage has been neglected. In course of this, the latter 
curves demonstrate the only zero bias anomaly discussed above, which is substantially 
smoothed, compared with the case of T = 4.2 K. On the contrary, the WKB scheme of 
calculation takes into account the variation of the potential profile inside the insulating 
barrier under applied bias voltage. In view of this, the differential conductances calculated 
by this scheme exhibit the tendency to increase at the direct bias voltage and to decrease 
at the reverse one. As it was shown above (see section 2.4), this behavior originates 
from the shift of the resonant levels inside the insulator due to externally applied electric 
field. This second mechanism in the origin of non-linear voltage dependence of impurity 
assisted conductance does not relate with the excitation of spin-flip processes. It becomes 
apparent at the voltages of the order of 50 mV and leads to the asymmetric voltage bias 
dependences in both cases of parallel and antiparallel configurations. 

Finally, the TMR amplitude as a function of the bias voltage is shown in Fig. 6 for 
the broad range of applied voltage. Its non-linear and asymmetric behavior in the range 
of 0.2 V originates from the asymmetry of the shifts of resonant impurity levels with 
respect to forward and inverse bias. The low bias voltage anomaly at 10 mV is also 
strongly pronounced at the curve corresponding to 4.2 K. The relative contribution of 
the impurity assisted conductance to the total current of electrons is diminished after the 
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value of applied bias voltage exceeds the half width of impurity band Ae/2 = 0.1 eV. 
Therefore, the TMR amplitude drops to value ~ 20% at 0.5 V corresponding primary to 
the pure tunnel conductance. 
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Appendix A 

Let w\ = c — a and W2 = b — c be the position of impurity with respect to the left and right 
interfaces, respectively, and w = b — a be the width of the insulator layer. We introduce 
tunneling densities of states for spin j (J,) electrons I^q) = fc^^o/l^l)^ + 9o) anc ^ 
denote 7| = |r-|-+|r_|_, 7^ = |r-f+|r^. Then, the posit ion- dependent quantum mechanical 
probabilities (c) can be found as follows, 
(a) In case of parallel configuration: 



TT(U) ( „\ _ MU) 



r 2 ul) e- 2 *> w /(w lW2 ) 

e -2q w 1 ~ j 

Ul)—^ — + r T(D" W2 



TT(U) /„\ _ ^ 1 TO) 



4 r 2 n e- 2 <^/( ra) 



M>2 

-TT(U) /-x _ 1 T 2 a) e~ 2 ^ w /( W1W2 ) 

-¥hr ' 9 e - 2 90»i e -2qow 2 ' l/ 1 ^ 

7l(T) W1 +7|(T)~ 

W2 



lu ,_ n , , _ 2 r T i>- 2 ^/( TO) 

7t U)i UIO 



o-_i(c) = <r_i(c) = - 



71- ^7L~ 



(b) In case of antiparallel configuration: 



r T r i e- 2go "'/(ii> im2 ) 

2(-D vw ~~ r e~ 2g0Wl , r e" 290 "' 2 
1 TU) - 5£ — + 1 KT) W2 



_TT(U) 

2 1 - 2> 



r T iy- 2 ^/(™) 



7TU)- 



-2(?oWi 



+ 7|(T)" 



-2g W2 



l«2 
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2^2' 



1 

9 



r r -2 qoW/{wiW2) 



-2gou>i 



7i(T)" 



wi 



+ 7TU)- 



u>2 



(A2) 



a nw)( z ) - 2 rf (i) e 2g0W /( W1W2 ) 



7t UIl 7j. ^2 



^UT)^) = ^ r T(l) e 290W /K^) 



Q -2g U>l -2q W2 ' 

71 ^7t~ 



The statistical probabilities Pj P {h) are independent of the configuration of the system. 
We denote h = pi B Hf/kT and Z = 2 cosh (2/i) + 1, then 



Pi 1 (fc) = pU(h) = Z~ x e lh , P[ 1 (h) = PH (/i) = Z~\ 

2 2 2 2 

P^T (fy = pU = Z -i e -2h > ( A3 ) 

2 2 



Appendix B 

The non-trivial functions i^, p are written as follows 
rIT a/ w eff, \ 



, (eV-2 / x fl gf)(e ey / fcr -l) 

1/21 ' 2 J ( e (ei/-2 MB fff )AT _ 1)(2 cosh(2/i B Ff /fcT) + 1) ' 

T U (v H e S) (eV + 2^ B Hf)(e^ kT -l)e^ kT 

^ ' 2 J ( e (eV+2 MB ^)/fcT -i)(2cosh(2// B #f/£;:r) + l)' 

,IT (eV + 2^f)(e^-l) 

-V 2 ^' 2 J (e( el/ + 2 ^^ ff )/ fcT -l)(2cosh(2/i B i/f/A;T) + l)' 

n s (eV-2» B Hf)(er v '« r -l)e-W"' 

' 2j (e(^- 2 ^^ cff )/ fcr -l)(2cosh(2/i B iJf/A;T) + l)' 



All the other ones, not written above, are equal to eV. 
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